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VOLUME IV: Ku~Band Antenna Acquisition and Tracking 

Performance Study 

VOLUME V: Hardware Simulation of Shuttle Ku-Band 

Antenna Acquisition 


This two volume report documents the work accomplished during 
amended Phase II of Contract Number NAS 9-14636. The system 
tradeoff analysis and system performance study that led to 
three antenna pointing and acquisition algorithms for hardware 
simulation is provided in Volume IV. The construction of 
results obtained from the hardware simulation are described 
in Volume V. A high level software description and a detailed 
software documentation of the computer programs is to be 
found in Volume V. The first three volumes of this report 
documents the work accomplished during Phase I of this contract 
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HARDWARE SIMULATION COMPUTER FLOW CHART 


OBJECTIVES 


LINCCMs overall objective under this contract was to develop 
a hardware simulation v;hich could emulate the Shuttle's Ku-Band 
Antenna Pointing and Signal Acquisition System. The desire was 
to develop a simulation in which the antenna pattern, the TORS 
search volume, the a priori probability distribution of TORS 
satellite position relative to the Shuttle, and the antenna scan 
procedure could be selected for the purpose of predicting 
performance for various Shuttle/TDRS antenna pointing and 
acquisition scenarios. 

The simulation was developed under the constaints of assuming 
a fixed network operation, the system must be real world imple^- 
mentable, it must be cost effective, the program execution time 
must not be excessive, the modulation technique is PN/BI-<f)/BPSK, 
and optimum performance is desirable. 




SUMMARY 


This document presents the results pertaining to the trade- 
off anaTysis and performance of the Ku-Band Shuttle antenna 
pointing and signal acquisition system. The study was perfoxmed 
assuming the existance of various antenna scanning, trajectories 
and various signal acquisition algorithms. The square, hexagonal 
and spiral trajectories were investigated assuming the TDRS 
postulated uncertainty region and a flexible statistical model' 
for the location of the TDRS within the uncertainty volume. 

The scanning trajectories, Shuttle/TDRS signal parameters and 
dynamics and three signal acquisition algorithms were integrated 
into a hardware simulation discussed herein and' documented in 
detail in Volume IV. The hardware simulation is- quite flexible 
in that it allows one to evaluate signal acquisition -performance 
for an arbitrary (programmable) antenna pattern, a large range 
of C/Nq*s, various TDRS/Shuttle a priori uncertainty distributions 
and three distinct signal search algorithms. 

Based upon the data made available during this contract 
period, certain Ku-Band forward link signal threshold characteristics 
were studied. The antenna pointing and acquisition, threshold Is- found 
to be less than 60 dB-Hz with an acquisition time dependent on the 
antenna scan procedure and acquisition algorithm implemented. 

Various techniques are discussed in this report and a computer 
program' is presented from which these can be evaluated. It 
appears that there will be no problem in meeting tlie system 
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spec even in light of the uncertainty associated v/ith the 

possibilities of antenna sidelobe acquisition. The details are 

provided herein. Assuming a single channel monopulse system and a 

30 MHz IF bandwidth, a monopulse tracking loop bandwidth of one Hz, 

the monopulse tracking jitter is 0.11 degrees at C/Ng = 54 dB-Hz. 

The Costas loop arm filter bandwidths can be chosen to be ten 

times the data rate such that no false-lock problem occurs 

during carrier acquisition. Assuming loop bandwidth of 

= 3 kHz, the loop jitter is 10 degrees at C/Ng = 60.4 dB-Hz. 

Acquisition can 'be accomplished in less, than 10 seconds. 

For the purpose of antenna scanning analysis, a reference 

coordinate system, whose z-axis is in line with the center 

axis of a specified uncertainty .cone of the TORS position, is 

chosen for the relative relation between the TORS and the Shuttle 

antenna. The scan path of the Shuttle antenna can be projected 

onto the (x,y) plane of the coordinate system; v/hile the uncertainty 

cone of the TORS position can be described by a circle in the 

(x,y) plane. The uncertainty in the position of the TORS is 

modeled by a truncated Gaussian probability density p(x,y) with 

2 2 2 . ' 

uncertainty parameters a = a = a . By changing the value 

x y 

of this variance parameter, the model is sufficiently general 
to include a uniform distribution of TORS position uncertainty 
((jp=«>) to one which specifies the position with probability one 
(yo). 

In the study of antenna scanning three types of trajectories 
are proposed and evaluated in terms of average scan time and the 
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structure of coverage over the uncertainty cone of the TORS 
position. The analysis technique used for finding average 
scan time is discussed and an illustration of the technique 
for a specific type of scanning trajectories and motion of the 
Shuttle antenna is given. The results show that the spiral tra- 
jectories is in general better than the other two trajectories, 
especially for a constant velocity along a trajectory, since the 
path of a spiral trajectory (from the center of uncertainty cone 
to its edge) can be shorter than the other two. 

From the Ku-band system specifications, it has been 
recognized that the variation of the received signal level at 
the Shuttle (due to TORS EIRP path loss and antenna pointing loss 
variations) varies as much as 23 dB. If the antenna sidelobes are 
not sufficiently suppressed, a potential problem called sidelobe 
acquisition may cause the degradation in system performance. 
Therefore in this tradeoff analysis study, three different 
acquisition strategies are proposed and evaluated. The first 
acquisition strategy may be used for the case where the sidelobe 
acquisition does not impose a problem to the system. The other 
two .strategies are primarily designed to avoid the sidelobe 
acquisition (especially the first sidelobe). 

In this study, we assume that the first sidelobe of the 
Shuttle antenna is 17.5 dB suppression from the main lobe of the 
antenna so that v/e can study the effectiveness of the acquisition 
algorithm proposed. Surely, one should note that the sidelobe 
acquisition can be overcome by tapering down the sidelobes of 
the Shuttle antenna below 23 dB. However, it has been 
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found that with one particular acquisition algorithm studied, 
a strong sidelobe of the antenna can speed spatial acquisition 
significantly. 

■The scan schemes and acquisition algorithms are all integrated 
by a computer simulation program. This software package, discussed 
■herein, provides a useful tool for predicting the performance 
of the KU'band antenna pointing system. Various options are 
available for users to do tradeoff studies on system parameters 
in designing an acquisition system. The software package has 
been tested on UNIVAC 1100 Series computers and verified with 
analytical results. 
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(1) Introduction 

The Ku-band antenna pointing acquisition systems based on 
the noncoherent signal energy detection and scanning scheme 
discussed previously have been simulated by digital computer^ 
for the tradeoff study in the design of the Shuttle antenna 
acquisition. Several acquisition algorithms implemented are 
given to provide flexibility in accommodating the different 
operating ranges of autotrack systems. The acquisition strategies 
are primarily oriented toward the avoidance of the first side 
lobe acquisition in the spatial search. In addition, the 
simulation also provides various options on the' scan tra'jector-ies, 
scan rates and a wide range of CNR. It essentially integrates 
the effect of interaction of spatial search acquisition schemes 
and the noncoherent signal energy detection. So the simulation 
program is able to support the design of hardware development. 

The simulation program was written in Fortran IV for the 
Univac 1100 series computer. Originally, it was .developed for 
time-shared operation, which is ’available at LINCOM Corporation. 
However, the program can also be used for batch processing with 
some modification. The program consists of 19 internal subroutine 
references and two external subroutines MERFIC (from IMSLS package) 
and BSSL (from MATHPACK) which are normally available for 
computer systems. 

The computing time to run the simulation program depends on 
the number of tests on the TORS position and the number of 
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samples taken for channel noises. These two variables are also 
provided' as options for a- user. The higher the values of these 
two variables, the higher the accuracy in predicting the acquisition •• 
performance and the longer the CPU time. The CPU time for a 
typical case is well within 10 sec if the number of noise samples 
is one and the number of simulated TORS position is less than 500. 

The simulation program provides, various options for users 
to test system performance for a wide range of system parameters. 

The following options included are; 

(1) The size of uncertainty cone of the TORS position and 
its statistical distribution function over the cone. 

(2) The antenna pattern of the Shuttle- receiver (3 dB 
beamwidth).'. 

(3) The types of antenna scan trajectories---square, hexa- 
gonal or spiral with different overlap,. 

(4) The scan rates--cohstant angular velocity, constant 
speed along a trajectory or combination of both’ previous 
cases (also a switching angle e^). 

(5) The probability of false-alarm and the operating range 
of carrier-to-noise ratio (CNR). 

(6) The receiver parameters — dwell time of Jthe-s.ignal energy 
detector, IF filter bandwidth, insertion loss and other 
circuit losses. 

■(7) The number of TORS positions used in the Monte Carlo 
simulation. 

(8) The numbev of samples taken for channel noises. 

(9) The types of acquisition algorithms. 
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(10) The choices of a deta'iled printout or final results 
in the simulation. 

The simulation program was structured as functional modules 
so that each subroutine can be easily replaced to adopt to 
any required function. For .instance, the antenna pattern used 
in the simulation can be simply replaced by other types of 
antenna pattern. To do so, one should remove the subroutine 
ANTNNA and substitute an equivalent subroutine ANTNNA. 

The simulation results are presented in two versions — one 
for a -time-shared mo.de and the other for a batch process mode. 

For the time-shared mode, the computer asks for input data 
through -interaction with users, A typical output generated by the 
simulation program consists of the following important parameters: 

(1) The threshold THl and TH2 used for the signal energy 
detector in a normal scan and in mini scans (such as 
adjacent and sidelobe scans), respectively, 

(2) Total number of cells to cover, the whole uncertainty 
region. 

(3) Maximum scanning time (= time spent to reach the last 
cell through normal path only). 

(4) Average scanning time and corresponding standard 
deviation. 

(5) "Probabilities" of detection, miss and false alarm, 

A sample printout is given in Fig. 1. 

(2) Functional Diagram', of Simulation Program 

To clarify the main functions performed in the simulation, 
a flow diagram of the simulation program i-s depicted in Fig, 2’. 
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RADIUS OF UNCERTAINTY REGION: 
BEAMW'IDTH: 

SPIRAL TRAJECTORY:: OVERLAP 

: CONSTANT ROTATIONAL SPEED: 
,^PROBABILITY OF FALSE ALA^RM: 

CNR EXPECTED RANGE: 

INTEGRATION TIME: 

IF FILTER BANDWIDTH: 

INSERTION LOSS: 

^OTHER LOSSES: 


10.0000 DEGREES 
1.6000 DEGREES 

30.000,0 % 

108.0000 DEGREES/S 
; 1000-02 • 

60.0000 DB ■ - 80.0000 DB 

5.0000 MS 

6. 0000 MHZ 
1.5000 DB 

.OOOOi DB- 


THE PROBABILITY OF FALSE ALARM WAS RAISED TO 0.1 TO AVOID SIDELOE 

.. .E DETECTION 
WHEN CNR = 80.00 DB/HZ 


THRESHOLDS: TH1 = 1.0568141 TH2 = 1.2276665 


NUMBER OF CELLS.: 302 


MAXIMUM SCANNING TIME = .3217689+02 SECONDS 


' NUMBER OF SATELLITE POSITIONS GENERATED: 
STANDARD DEVIATION: 


500 

. 1000000+01 


CARRIER-TO-NOISE ' RATIO: 

ACQUISITION AVERAGE TIME: 

STANDARD DEVIATION: 

PROBABILITY OF DETECTION: 

PROBABILITY OF DETECTION (SL. VAR.;: 

PROBABILITY OF HISS: 

PROBABILITY OF FALSE ALARM (IN THE SIMULATION): 


..'6'000000+02 DB 
.3023562+0,1 S 
,.2362385+01 S 
. 1000000+01 
.0000000 
.-0000000 
.0000000 


Figure Ta. A Sample of Input and Output Data. 
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The TORS position is simulated by a random generator with a. 
specified distribution function. To measure the spatial 
acquisition time of the TORS position, a set of s 'Emulated TORS 
positions is tested for a specified trajectory and scan rates 
(speed for constant angular velocity and/or constant velocity 
along a trajectory). The simulated data reported here are .all 
obtained by 500 simulated TORS positions,. 

In the simulation program, the acquisition monitor is not 
a separate unit, but imbedded in the subroutine ACQ. This is 
implemented just for the convenience in programming. The 
detailed structure of the program is referred’ to -herein. 

(3) Capabilities of Software Package 

The program developed for the Ku-band antenna pointing system 
provides- various options in the size of uncertainty cone, for the 
TDRS positions, the type of scanning and trajectories-,, and spatial 
acquisition algorithms. Here are the capabilities of the software 
package: 


(1) 

( 2 ) 

(3) 

■(4) 

(5) 


( 6 ) 


Conical angle of the TDRS position can be 4°, 8° or 10°. 
2 

The variance of a priori probability distribution 
function of the TDRS position can vary from 0.1° to 10°.' 
Antenna pattern can have beamwidth ^ 10°. ' 

Three scan trajectories and wide range of' scan rates 
(but limited by electrical dwell time) are available. 

The noncoherent signal detector can be operated over 
a wide range of CNR and Bj^T .-product greater than 1000, 
Spatial acquisition algorithms with/wi^-hout avoidance 
of si delobe acquisition are selectable-. 
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(7) Probability of false alarm can be from 0.5 down to 10~^. 

(8) The number of TORS positions in Monte Carlo simulation 
can be selected by users. But it is limited by' the CPU 
time in executing the simulation program. 

• In addition, the software package also provides a useful 
feature for system designers to interact with a computer if the- 
time-shared facility is available. Upon user's choices, a detailed 
spatial acquisition procedure can be printed out for a closed 
examination on the behavior of acquisition systems. Finally 
one should also note that the application of the software 
package should not be limited to the Ku-band system parameters 
even though it is oriented tov/ard the study of Ku antenna pointing 
systems. 

(4) Computer Program Utilization 

The simulation program was- designed for two modes of usage — 
batch process and time-shared modes. Since the computer program 
was originally developed through time-shared mode, it provides the 
mechanism of interactions between computers and users. However, this 
set of interactions should be suppressed for the batch processing. 

In order to run the program in batch, one should provide the same set 
of input data as one .type in input data in the time-shared mode’, 
except for the first input data. To illustrate the details of the 
input data required, several examples are given in Figs. 1 and 2 
to show the options available in the simulation and the corresponding 
sets of input data for the batch process are also listed. 

(1) Example 1 

(a) Time Shared Mode ' 
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IF YOU ARE OiJ-LlNE, PLEASE TYPE 1 
>1 

*** ON-LINE *«< 




INPUT DATA (F.REE FORMAT) 

(REAL) RADIUS OF UNCERTAINTY REGION (DEGREES) 

(REAL) EEAHWIDTH (DEGREES) 

(INTG) TYPE Or TRAJECTORY: KSQUARE) 2(HEXAGONAL) 3(SPIRAL) 

> 10 . , 1 . 6,1 ' 

(INTO) TYPE OF MOTION: 

1 - CONSTANT ANGULAR VELOCITY 

2 - CONSTAfJT VELOCITY ALONG TRAJECTORY 

3 - COMBINATION OF "1" AND "2" 

>2 

(REAL) ANGULAR VELOCITY (DEGREES/S) 

(REAL) VELOCITY ALONG TRAJECTORY (DEGREES/S) 

>106. ,5.0Mi5 

(REAL) PROBABILITY OF FALSE ALARM 

(REAL) EXPECTED MINIMUM CARRIER-TC-NOISE RATIO (DB) 

(REAL) EXPECTED MAXIMUM CARRIER-TO-KOISE RATIO (DB) 

>i.e-ii,eo.,eo. 

CNR VALUES TO BE USED IN THE SIKULATICfi 
(REAL) (MINIMUM, MAXIMUM, INCREMENT. GE . 0. ) 

>60. , SO. ,5. 

RECEIVER PARAMETERS: 

(REAL) INTEGRATION TIME (MS) 

(REAL) IF FILTER BANDWIDTH (MHZ) 

(REAL) INSERTION LOSS (DB) 

(REAL) OTHER LOSSES (DB) 

>5. , 6. , 1 .5, 0. 

(INTG) NUMBER OF SATELLITE POSITIONS TO BE GENERATED DURING THE SIMULATION 
>500 

(INTG) 0 - SATELLITE POSITIONS ARE ENTERED AS INPUT DATA 
1 - SATELLITE POSITIONS ARE INTERNALLY GENERATED 

j,., . - — — - 

(REAL) STANDARD DEVIATION VALUES: (MINIMUM, MAXIMUM, IHCRE-MENT. GE . 0. ) IN DEGREES 

>2. 5, 2. 5,1. 

(INTG) INTERMEDIATE RESULTS: NC(0) YESd) 

>0 


TRANSITION MODES AND VERIFICATION TIMES: 

(REAL) NORMAL TO SIDELOBE SCANNING 
(REAL) NORMAL TO SURROUNDING CELLS SCANNING 
(REAL) VERIFICATION TIME 
>1.e-3,1.e-3,5.c-3 

SCANNING SCHEME: 

(INTG) .0 - N0RMA1.--_P.ATH ONLY 

1 - NORMAL PATH AND SIDELOBE CELLS 

2 - NORMAL PATH, NEARBY, AND THEN SIDELOBE CELLS 

>0 

(INTG) NUMBER OF NOISE SAMPLES (EVEN K(.)<20 ) 


Figure 3. A Sample Interaction Between a Computer and User. 
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IF YOU ARE OK-LIKE, PLEASE TYPE 1 
>1 

*'*« OK-LIUE *»* 
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INPUT DATA (FREE rCRKAT) 

(REAL) RADIUS OF UNCERTAINTY REGION (DEGREES) 

(REAL) EEAEWIDTh’ (DEGREES) 

(IKTG) TYPE Or TRAJECTORY: KSQUARE) 2{HEXAGOHAL) 3(SPIRAL) 

>10. .1.6,3 

(REAL) OVERLAP ( %-) 

>30. ■ 

(INTO) TYPE OF KOTICN; 

1 - CONSTAI.T ANGULAR VELOCITY 

2 - CONSTANT VELCCllY ALONG TRAJECTORY 

3 - COMBINATION OF "I" AK=P "?« 

>3 

(REAL) ANGULAR VELOCITY (DLGREES/S) 

(REAL) VELOCITY ALONG TRAJECTORY (DEGHEES/S) 

(REAL) SWITCHING ANGLE (DEGREES) 

>108. ,5.0NiL,5.5 


(REAL)' PROBABILITY Or FALSE ALARM 

(REAL) EXPECTED MINIMUM CARRIER-TO-NOISE RATIO (DB) 

(REAL) EXPECTED MAXIMUM CARRIER-TC-NOISL RATIO (DB) 

►1.e-3,60.,80. 

CNR VALUES TO EE USED IN THE SIMULATION 
(REAL) (MINIMUM, MAXIMUM, INCREMENT. GE. 0. ) 

*75.,75.,1. 

RECEIVER PARAMETERS: 

(REAL) INTEGRATION TIKE (MS) 

(REAL) IF FILTER BANDWIDTH (MHZ) 

(REAL) INSERTION LOSS (DB) 

(REAL) OTHER LOSSES (DB) 

’5. , 6. , 0.5 1 1 . 

(IHTG) NUMBER OF SATELLITE POSITIONS TO EE GENERATED DURING THE SIMULATION 
■1000 

(INTG) 6 - SATELLITE POSITIONS ARE ENTERED AS INPUT DATA 
1 - SATELLITE POSITIONS ARE 7»TrR^’ aMY GPNPP-CrJi: 

•1 

(REAL) STANDARD DEVIATION VALUES: (MINIMUM, MAXIMUM, INCREMENT. GE. 0. ) IN DEGREES 

>1.0, 5. 5, 1.5 

(INTO INTERMEDIATE RESULTS: N0(0) YESd) 

■1 

(INTG) HOW MANY RESULTS DO YOU WANT TO SEE IN DETAILS? 

■30 


TRANSITION MODES AND VERIFICATION TIMES: 
(REAL) NORMAL TO SIDELOBE SCAIH-’ING 
(REAL) NORMAL TO SURROUNDING CELLS SCANNING 
(REAL) VERIFICATION TIME 
■1 . e-3 » 1 • c-3, 5 .e-3 


SCANNING SCHEME: 

(INTG) 0 - NORMAL PATH ONLY 

1 - NORMAL PATH AND SIDELOBE CELLS 

2 - NORMAL PATH. NEARBY. AND THEN SIDELOBE CELLS 
■2 ' , 

(INTO NUMBER OF NOISE SAMPLES (EVEN K(.)<20 ) 
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IF YOU AHE ON-LINE, PLEASE TYPE 1 

>1 

*** ON-LINE 
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INPUT DATA (FREE FORMAT) 

(REAL) RADIUS OF UNCERTAINTY REGION (DECREES) 

(REAL) EEAKWIDTH (DEGREES) 

(INTO TYPE OF TRAJECTORY: 1 (SQUARE) 2 (HEXAGONAL) 3 (SPIRAL) 

> 10 . , 1 . 6, 2 ' ■ 

(INTO TYPE OF MOTION: 

1 - COr.'.S~flHT funu-.a^ VfinrT.IY, 

2 - CONSTANT VELOCITY ALONG TRAJECTORY 

3 - COMBINATION OF "I" AND "2" 

>1 

(REAL) ANGULAR VELOCITY (DEGREES/S) 

>108. 

(REAL) PROBABILITY OF FALSE ALARM 

(REAL) EXPECTED MINIMUM CARRIEF-TO-NOISE RATIO (DB) 

(REAL) EXPECTED MAXIMUM CARRIER-TO-NOISE RATIO (DB) 

>1.e-iJ,60.,80, 

CNR VALUES TO EE USED IK THE SIMULATION 
(REAL) (KTKIMUM, MAXIMUM, INCREMENT. GE . 0. ) 

>60. ,80, ,5. 

RECEIVER PARAMETERS: 

(REAL) INTEGRATION TIME (MS) 

(REAL) IF FILTER BANDWIDTH (MHZ) 

(REAL) INSERTION LOSS (DB) 

(REAL) OTHER LOSSES (DB) 

> 5 . , 6 . . 1 . 5 , 0 . 

(INTO NUMBER OF SATELLITE POSITIONS TO BE GENERATED DURING THE SIMULATION 

>10 

(INTO 0 - SATFT.1.TTE PORITTAMc; tpc FHTERED AS INPUT DATA ' 

1 - SATELLITE POSITIONS ARE INTERNALLY GENERAIED • 

>0 

(INTO INTERMEDIATE RESULTS: N0(0) YES(1) 

>0 


ORIGINAL PAGE IS 
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(INTO NUMBER OF NOISE SAMPLES (EVEN K(.)<20 ) 

>N 


TRANSITION MODES AND VERIFICATION TIMES': 

(REAL) NORMAL TO SIDELOBE SCANNING 
(BEAL) NORMAL TO SURROUNDING CELLS SCANNING 
(REAL) VERIFICATION TIME 
>1.e-3,1.c-3,5.e-3 

SCANNING SCHEME: 

(INTO) 0 - NORMAL PATH ONLY 

1 - NORMAT. PATH AKD SIDELOBE CELLS 

2 - NORMAL MTH, KEAHbTT^ND THEN SIDELOBE CELLS 

>1 



Figure 3 , (Cont'd-) 
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(5) Functional Descriptions o'f Subroutines 

The functions of each subroutine used in the computer program 
for UNIVAC 1100 series are briefly described one-by-one here. 

Their flow charts and" detailed documentations are given in 
Volume II separately. The description of subroutines is 
given in alphabetical order to be consistant with the compu>ter 
printouts from UNIVAC 1.108 time-shared mode. The program is, 
in ict, started from MAIN then SIMULA. The subroutine SIMULA, 
the core of the simulation, cells andmonitors the operations 
of the Monte Carlo simulation- 
Subroutine. ACQ 

This s.imulates the acquisition algorithms discussed 
previously. The variable NULL, used to select the type^of 
the acquisition algorithms for the simulation, is defined, as 
follows: 

NULL = 0 for normal scan 

1 for normal scan plus sidelobe scan 

2 for normal scan, adjacent scan and 
sidelobe scan 

In the process of the spatial acquisition, the time spent in 
any scan, including the transition from one scan to another 
and verification time, is accumulated as one of the outputs of 
the subroutine. The other important output parameters of the 
subroutine are: 

(1) The status of the spatial search — success or failure 
(IPMISS and IPDET). 

(2) How to acquire--the detection- was declared by.- a normal 
scan or pi us sidelobe scans (LPDSL). 
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(3) The acquired location of the TORS (KHIT (IHIT)-). 
Subroutine ANTNNA 

The antenna, pattern of the simulated Shuttle receivers is 
assumed to have a first sidelobe of -175 dB from the peak of the 
main lobe. The antenna gain is computed from the offset angle 
between the simulated TORS position and the- boresi.ght axis of 
the Shuttle. Hence, the offset angle (TH) is an input variable 
in addition to the parameters characterizing the antenna (THB, 
PATTC). A control variable K is used for computing the 
antenna parameters as shown here 


K 


-2 for computing the angle of first null 

-1 for computing the angle of the peak of 

of the second sidelobe 

0 for computing the angle .of the peak of the 
first sidelobe 

+1 for -3 dB beamwidth 

2 for normal usage in computing the gain 
of the antenna at offset TH. 


Subroutine CELL 

For a specified trajectory (KSCAN), the subroutine computes 
the locations of search cells (XCTR,YCTR) along the trajectory 
and the scan time to reach each search cell from the designated 
center of the uncertainty zone of the TORS position. The scan 
time is computed by calling the subroutine TCENTR. Hence, it 
requires the following input variables 


KSCAN = 


POOR QTJALTIi' 


V . 


1 for square trajectory 

2 for hexagonal trajectory 

3 for spiral trajectory 
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MOTION = 1 for con'stant rotational speed 

2 for constant speed along trajectory 

3 for constant rotational speed first 
then constant speed along trajectory 

VANG = angular velocity for constant rotational 
speed (in radian/sec) 

VLIN = constant speed along trajectory (in radian/sec) 

To ensure the electrical dwell time to be smaller than times 
of the minimal mechanical dwell time, a warning is printed out 
if the antenna scan is too fast. In the subroutine is set 
to be.ljj 

Subroutine GAUSS 

This is a Gauss random generator whose mean and standard 
deviation are specified by AVR and SIGMA, respectively. The 
Gaussian samples are obtained in pair by the following simple 
relations 

1 /2 

y^ = a(-2.1n x-j ) ‘ cos 2 itX 2 + u 

1 /2 

y 2 .= cf(-2.1n x^) ^ sin 2 uX 2 + y 

where x-j and X 2 a pain of independent random variables uniformly dis- 
tributed between 0 and 1 , and y and o are its mean and variance, 
respectively. -The algorithm has been tested for its mean, variance 
and skew coefficient. The results show that the al gorithm gives 
a high confidence in its statistical nature. 

Subroutine HIT 

The signal energy detector is simulated by this subroutine. 

The received signal energy (sum of signal power and noise power) 
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is compared 'wi th a preset threshold (Thi). If the threshold is 
crossed, a hit is declared (HIT=.TRUE. }; otherwise, .no hit is given. 

The input variables are the position of the TORS and the location 
of the boresight axis of the Shuttle antenna. The offset angle 
TETA is computed, then the antenna gain by calling subroutine 
ANTNNA, the signal power and noise power. 

Subroutine HPOSTN 

This subroutine computes the centers of hexagonal cells along 
the hexagonal trajectory. The computation is proceeded from 
cell to (K+1)^ cell. 

Subroutine ITR 

This subroutine implements the iterative algorithm, needed 
in solving for eqs. (20) and (21). The control variable K 
is used to indicate the application of the algorithm; to- eq. (20) 
or eq. (21). 

MAIN 

The function of the MAIN defines the dimension of variables 
used in' the simulation. The program has been set to deal with 
the maximal number of search cells to be 500.. If the number- 
of the search cells goes beyond 500, one should expand dimensions 
of all variables in the first dimension statement', 

MAIN/HAP 

This is a set-up for executing the simulation program. Two 
system subroutine packages — MATHPACK and . IMSL-rare called to 
satisfy the external references. 

Subroutine NOISE 

The simulated channel noise power at the -output of the- integrator- 







is computed here. The number of samples for equivalent Gaussian 
noises is set , by the variable KNOISE. 

Subroutine QPOSTN 

It computes the center of (K+1)— search cell from the known 
"t h 

center of the K— search cel l along, the square trajectory.. 
Subroutine SIMULA 

This is the core of the simulation program. It reads in the 
■data,' prints the outputs and also monitors the operation of the 
program. The detailed operation of the subroutine is referred to 
herein. 

_Subroutine 5LCELL 

The search cells alons a sidelobe trajectory and adjacent 
cells to a hit cell are computed here. The adjacent cells are 
structured as those for hexagonal cells. Hence-j the 'number of 
the adjacent cells to be scanned is fixed as six. However, 
the number of cells for the .sidelobe trajectory varies depending 
upon the type of trajectory used for normal scan. 

Subroutine SP05TN 

* s t 

It computes the center of (K-H) — cell along a spiral 

trajectory from the known center of the K— cell . 

Subroutine TCENTR 

It computes the scan time required' for the Shuttle antenna 
to sweep along a trajectory to reach a spatial point in ^ the 
uncertainty zone of the TORS position. It requires to input 
the type of motion, scan rate (constant rotation and/or constant 
speed along a trajectory) and the coordinates o-; the point. 


ORIGINAL PAGE IS 
OP POOR QUAUTZ 
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Subroutine TORS 

The subroutine generates a simulated TORS position according 
to a specified Gaussian distribution. If a Gaussian sample falls 
outside the uncertainty zone of the TORS position, it is discarded 
and another sample is taken. Its input variables are SIGMA, 

THM2, and its output variables XTDRS and YTDRS. 

Subroutine TH1TH2 

The setup of threshold levels used for spatial -acquisition 
is done in this subroutine. The threshold THl is first computed 
for the given probability of false-alarm under the channel noise 
condition. Then it 'checks the- probability of detection 

T3.S 

due to the second sidelobe of the Shuttle antenna at the 
maximal value of CNR in the specified operating range. If 
is larger than the given probability of false alarm, the threshold 
is adjusted to yield at least 0.1 probability of false alarm 
due to second sidelobe at the strongest CNR. This implies that 
the probability of false-alarm at lowest CNR is much smaller 
than the specified one. 

The threshold TH2 is computed based on the assurance of 
having the specified probability of false alarm due to the first 
sidelobe at strongest CNR. 

Subroutine UNIFOR 

This is a uniform random generator over (0,1). A pair 
of output samples are -generated for each cell. 
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APPENDIX I 

SUBROUTINE DESCRIPTIONS 


1 . MAIN -Program 

Purpose: Define the working areas needed throughout 

the simulation. 


2. Subroutine SIMULA 


Purpose: Read input data, determine the "a priori" 

set of data necessary for the "acquisition 
process," and compute desired time averages 
and probabilities. 


Usage: Call SIMULA{NCELL ,XCTR,YCTR,TCTR,ELEV'',AZIM,NSLDIM, 

XSLO , YSLO , XSL ,YSL ,NNEARD , XNEARO , YNEARO , 
XNEAR,YNEAR,KHIT,TH) 

INPUT PARAMETERS 

NCELLD - MAXIMUM NUMBER OF CELLS .POSSIBLE TO' BE- USED TO SC-AN THE 

UNCERTAINTY REGION 

HSLDIH - MAXIMUM NUMBER OF CELLS POSSIBLE TO BE USED IN THE 

SIDELOBE SCANNING 

NNEARD - MAXIMUM NUMBER OF CELLS POSSIBLE, TO BE USED IN THE 

"SUHROUHDING CELLS" SCANNING 

OUTPUT PARAMETERS 

(XCTR,YCTR) - VECTORS WITH THE CENTERS OF THE VARIOUS CELLS 
(CARTESIAN COORDINATES) 

• (ELEV,AZIM) - VECTORS WITH THE -CENTERS OF THE VARIOUS CELLS 

(POLAR COORDINATES) 

(XSLO, YSLO) - VECTORS VMTH THE CENTERS OF THE SIDELOBE CELLS WHEN 
ANTENNA BORESIGHT IS POINTING TO THE ORIGIN 

-(XNEARO, YNEARO) VECTORS WITH THE CENTERS 'OF THE NEARBY CELLS WHEN 
ANTENNA BORESIGHT IS POINTING TO THE ORIGIN 

KBIT - VECTOR WITH THE LAST HITS, STARTING WITH THE LAST CELL 

IN THE NORMAL PATH WHERE A HIT WAS OBSERVED 

TH - VECTOR WITH THE THRESHOLDS USED .TO - SEARCH NEMBY CELLS 
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INTERNAL VARIABLES 
THH 
THB 
OVLAP 
KSCAN 


MOTION 


VANG 

VLIN 

SWTCH 

PFA 

CNRMIN 
- CNRHAX 
CNR 


RADIUS OF THE UNCERTAINTT REGION (RADIANS) 

RADIUS OF -3DB ANTENNA MAIN LOBE- 
OVERLAP 

TYPE OF TRAJECTORY 

1 - SQUARE 

2 - HEXAGONAL 

3 - SPIRAL 

DESIRED ANTENNA MOTION . 

1 - CONSTANT ROTATIONAL SPEED 

2 - CONSTANT VELOCITY ALONG TRAJECTORY 

3 - COMBINATION OF "1" AND ”2" 

ANGULAR VELOCITY (RADIANS/S) 

LINEAR VELOCITY ALONG TRAJECTORY (RADIANS/S) 
SWITCHING ANGLE (IN RADIANS) TO GO FROM "MOTION=V* 
TO "MOTION =2" 

PROBABILITY’ OF FALSE ALARM 

EXPECTED CARRIER-TO-NOISE RATIO DYNAMIC RANGE (DB) 
ACTUAL CARRIER-TO-NOISE RATIO (IN DB) USED IN THE 
SIMULATION 


TD 

BIF 

FLOSS 

OTHERL 

NTDRSS 

SIGHAP 


TSL 

TCLOSE 

TVERIF 

NULL 


TH1.TH2 

PDET 

PDSL 

PMISS 

PRPA 

TAVR 

SIGHAT 


ELECTRICAL DWELL TIME (MS) 

- IF FILTER BANDWIDTH (r:H2) 

FILTER INSERTION LOSS (DB) 

OTHER LOSSES (DB) 

NUMBER OF SATELLITE POSITIONS GENERATED DURING THE 
SIMULATION FOR EACH "SET OF CONDITIONS" 

STANDARD DEVIATION (IN RADIANS) OF THE TRUNCATED 
GAUSSIAN -DENSITY FUNCTION USED 'TO GENERATE THE SATELLITE 
POSITION 

TIME INTERVAL TO SWITCH FROM NORMAL SCANNING TO SIDELOBE 
SCANNING 

TIME INTERVAL TO SWITCH FROM NORMAL SCANNING TO NEARBY 
SCANNING 

VERIFICATION TIME 
SELECT SCANNING SCHEME 

0 - NORMAL PATH ONLY (AUTO-TRACK RANGE < ANGLE BETWEEN’ 

1 - NORMAL PATH AND SIDELOBE CELLS (AUTO-TRACK RANGE 

< ANTENNA NULLS 

2 - NORMAL PATH, SIDELOBE, AND THEN NEARBY CELLS (AUTO- ‘ . 

TRACK RANGE < BEAMWIDTH) 

3 - NORMAL PATH, NEARBY, AND THEN SIDELOBE CELLS (AUTO- 

TRACK RANGE < BEAMWIDTH) 

THRESHOLDS 

PROBABILITY OF DETECTION 

PROBABILITY OF DETECTION BY SIDELOBE ACQUISITION 
PROBABILITY OF MISS 

PROBABILITY. OF FALSE ALA'.RM (SIMULATED) 

AVERAGE TIME (IN SECONDS) TO DETECT THE SATELLITE 
CORRESPONDING STA'NDARD DEVIATION 


OEIGINAL PAGE JS 
OP POOR QUAury 





Figure'Al. Flow Diagram for Subroutine SIMULA. 
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DETERMINE 

T - TIME SPENT TO DETECT SATELLITE 
(XDET,YDET)-FINAL POSITION OF SHUTTLE 
ANTENNA BORESIGHT. 

IPDET JPDSLJPMISS 


1 

TAVR=TAVR+T 

SIGMAT=SIGMAT+T*T 



Fig. A1 . Continued 
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Fig. A1 . Continued 
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Fig. Al. Continued 
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16. Subroutine HIT 


Purpose:- Simulate the noisy output of the integrator and 

verify if this output is greater than a prescribed 

threshold (HIT=TRUE) OR NOT (HIT=FALSE). 

Usage: Call HIT(X,Y,XTDRS,YTDRS,TH,RECSIG,HIT,S0) 

Input Parameters': 

(X,Y) -ANTENNABORESIGHT- 

(XTDRS,YTDRS)~ SATELLITE POSITION 
TH - THRESHOLD 

Output Parameters: 

RECSIG - RECEIVED SIGNAL POWER + NOISE 

HIT - LOGICAL VARIABLE 

TRUE - IF RECSIG > TH 
FALSE - IF RECSIG < TH 

SO - VARIABLE USED IN THE NOISE GENERATION 

Subroutines Required: NOISE 

17. Subroutine NOISE 


Purpose: Generate the equivalent noise power at the output 

of the integrator. 

Usage: Call NOISE{S0} 

Input Parameter: 

SO - VARIABLE USED IN THE GAUSSIAN NOISE GENERATION 

SIGMAN - GAUSSIAN NOISE STANDARD DEVIATION (COMMON) 

POTXG - RECEIVED SIGNAL -POWER (COMMON) 

-NSAMPL - IF FILTER BANDWI-DTH*INT£GRATION TIME (COMMON) 

Output Parameter: 

EQN - EQUIVALENT NOISE AT THE OUTPUT OF THE INTEGRATOR 

(COMMON) 

Subroutines Required: GAUSS 

Model : 

ORIGINAL PAGE IS' 

SIGN2 = SIGMAN*SIGMAN OP POOR OUAT TW 

AVR = SIGN2 yUAUiY 

STD = SQRT((l.+2.*POTXG/SIGN2)/NSAMPL)’*SIGN2 
EQN = N( AVR, STD) 
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14. Subroutine UNI FOR 


Purpose: Generate two uniformly distributed random numbers 

•in the interval (0,1). 

Usage: Call UNIF0R(S,X1 ,X2) 

Input Variable: 

S - INPUT RANDOM NUMBER (MODIFIED IN THE SUBROUTINE) 
Output Variables: 

XI ,X2 - PAIR OF INDEPENDENT RANDOM VARIABLES UNIFORMLY 

DISTRIBUTED BETWEEN 0 AND 1 


Method : 

X =M0D(S*7**5, 2**21 -1)/2**31 
-15. Subroutine ACQ 


Purpose: Simulate the Ku-Band acquisition system. 

Usage: 

Call 


ACC (NULL, TUB, XKS, NCE.LL, XCTR ,YCTR , TCTH , AZIM , MSL, 

XSLO.YSLO, DTSL, NNEAR, XNEARO, YNEARO.DTN'EAR, VANG, 
TH1 ,TH2,TSL, TCLOSE,TVERIF ,XTDRS, YTDHS, NNEARD, SO, 
PRT,T,IPKISS, IPDET, IPDSL.XDET, YDET.KHIT, IHIT, TH, 
XSL,YSL,XNEAR,YNEAR) 
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INPUT PARAMETERS: 
NULL 


THB 

XKS 

NCELL 

(XCTR.YCTR) - 

TCTR 

A2IM 

NSL 

(XSLO.YSLO) - 
DTSL 


SELECTS SCANNING SCHEME 

0 - SEARCH CELLS IN THE NORMAL PATH ONLY WITH 

THRESHOLD Till 

1 - SEARCH CELLS IN THE NORMAL PATH WITH THRESHOLD- TH1 

AND THEN SIDELOBE CEJ.LS WITH THRESHOLD TH2 
(AUTO-TRACK RANGE < ANGLE BETWEEN NULLS) 

3 - SEARCH CELLS IN THE NORMAL PATH WITH THRESHOLD TH1 , 
THEN NEARBY CELLS WITH THRESHOLD, AND THEN SIDELOBE 
CELLS WITH SAME THRESHOLD TH2 
(AUTO-TRACK RANGE < ANGLE BETWEEN NULLS) 

RADIUS OF ANTENNA MAIN LOBE 
DISTANCE BETWEEN SUCESSIVE CENTERS 

NUMBER OF CELLS TO SCAN THE WHOLE UNCERTAINTY REGION 
VECTORS WITH THE CENTERS OF THE CELLS 

VECTOR WITH THE TIME INSTANTS EACH CENTER IS ACHIEVED 
VECTOR WITH THE AZIMUTH OF EACH CENTER 
NUMBER OF SIDELOBE CELLS 

VECTORS WITH THE CENTERS OF THE SIDELOBE CELLS 
WHEN THE ANTENNA IS POINTING TO THE ORIGIN 
TIME INTERVAL TO GO FROM ONE SIDELOBE CELL TO THE 


NNEAR 

(XNEAROjYNEARD) 

DTKEAR 

VANG 

TH1,TH2 

TSL 

TCLOSE 


NEXT ONE 

NUMBER OF NEARBY CELLS 

VECTORS WITH THE CENTERS OF THE NEARBY CELLS 
TIME INTERVAL TO GO FROM ONE NEARBY CELL TO THE NEXT 
CONSTANT ROTATIONAL SPEED 
THRESHOLDS 

TIME INTERVAL TO GO FROM 
THE SIDELOBE VARIATION 
TIME INTERVAL TO GO FROM 


ONE 


A NORMAL SCANNING SCHEME TO 
A NORMAL SCANNING SCHEME TO THE 


TVERIF 

(XTDRSjYTDRS)- 

NNEARD 

SO 

PRT 


NEARBY VARIATION 
VERIFICATION TIME 
GENERATED SATELLITE POSITION 

DIMENSION OF VECTORS XNEARO, YNEARO, XNEAR , YNEAR, 
AND TH IN THE CALLING PROGRAM 
VARIABLE USED FOR RANDOM NOISE GENERATION 
LOGICAL VARIABLE 

TRUE IF A PRINTOUT SAMPLE IS DESIRED 


OUTPUT PARAMETERS 
T 

IPMISS 
IPDET 
- IPD.SL 

KHIT 

IHIT 

TH 


TIME SPENT TO DETECT THE SATELLITE 

(0) THE SATELLITE WAS DETECTED 

(1) THE SATELLITE WAS MISSED 

(0) THE SATELLITE WAS HISSED 

(1) THE SATELLITE WAS DETECTED 

(0) THE SATELLITE WAS NOT DETCTED BY SIDELOBE DETECTION 

(1) THE SATELLITE WAS DETECTED WHEN THE SYSTEM WAS 
SEARCHING SIDELOBE CELLS ' 

VECTOR WITH THE LAST HITS, STARTING WITH THE LAST CELL 
IN THE NORMAL PATH WHERE A HIT WAS- OBSERVED 
KHIT(IHIT) DEFINES THE CELL WHERE THE SATELLITE WAS 
DETECTED 

VECTOR WITH THE THRESHOLDS USED TO SEARCH THE NEARBY 
CELLS CORRESPONDING TO THE LAST CELL SEARCHED IN THE 
NORMAL PATH 


INTERNAL VARIABLES: 

(XSLjYSL) - VECTORS WITH THE CENTERS OF SIDELOBE CELLS 

CORRESPONDING TO A GENERIC POSITION IN THE NORMAL PATH 
. (XNEAR, YNEAR) - VECTORS WITH THE CENTERS OF NEARBY CELLS 

CORRESPONDING TO A GENERIC POSITION IN THE NORMAL PATH 


ORIGINAli 

OP POOR QI I ATJa? ..... 
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Subroutines Required: 

ANTENNA, TH1TH2,SLCELL,UNIF0R,TDRSS,ACQ. 
Function Di-agram (see Fig. A2) 
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IP = 2 
TH(IP)=TH2 
T=T+TVERIF 


HIT-:TH2 



HIT: TH2 



OEIGINAt PAGE IS 

OE PooB quality; 


Fig. A2. Continued 
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Figure A2. Continued. 
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3. Subroutine ANTNNA 


Purpose: Determine angle and gain corresponding to ' 

(a) first null, (b) second sidelobe, (c) first 
sidelobe, (d) -3 dB point in the mainlobe, and 
(e) gain for any input angle.. 

Usage; Call ANTNNA (K,TH,GAIN) 

Input Parameters: 

K - CONTROL PARAMETER 

-2 - THE SUBROUTINE DETERMINES ANGLE FOR FIRST ZERO 
-1 - THE SUBROUTINE DETERMINES ANGLE FOR SECOND MAXTMUM' 

AND CORRESPONDING GAIN 

0 - THE SUBROUTINE DETERMINES' ANGLE FOR FIRST MAXIMUM 

AND CORRESPONDING GAIN 

1 - THE SUBROUTINE DETERMINES ANGLE CORRESPONDING ' 

TO MAINLOBE -3dB POINT AND CORRESPONDING GAIN 

2 - THE SUBROUTINE DETERMINES THE GAIN FOR THE 

ANGLE = TH 

TH - INPUT ANGLE (K=2) ' 

Output Parameters: 

TH - ANGLE OF FIRST NULL (K=-2), OR SECOND SIDELOBE 
• (K=-T), OR FIRST SIDELOBE (K=0), OR MAINLOBE 
' -3 dB POINT (K=l) 

GAIN - ANTENNA GAIN 


Subroutines Required: 

Bessel 1 Function: BSSL (MATHPACK) 

Model : 

GAIN = J-|(x)/x; X = sin'^(TH) + 1.61/THB 
4. Subroutine TH1TH2 

Purpose: Determine thresholds THl and TH2 for a given 

probability of false alarm. 

Usage: CALL THl TH2 (CN21 ,CN22,PFA,BSF, FLOSS, OTHERL, 

NSAMPL,TH1,TH2) 


Input Parameters: 
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CNRl 

PFA 

BIF 

FLOSS 

ORTHERL 

NSAMPL 


CNR2 - CARRIER-TO-NOISE DYANMI.C RANGE (DB) 
-PROBABILITY OF FALSE ALARM 
-IF FILTER BANDWIDTH 
-FILTER INSERTION LOSS (DB) 

-OTHER LOSSES (DB)’ 

-NUMBER OF SAMPLES 


Output Parameters: 

TH1 -NORMAL PATH SCANNING THRESHOLD 

TH2 -SIDELOBE CELLS SCANNING THRESHOLD 

Subroutines Required: 

ANTNNA,MERFCI (IMSL) 

Function Diagram (see Fig. A3) 


4T 
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5. Subroutine CELL 

Purpose: Compute number of cells needed for complete 

scan of the uncertainty region, their centers 
and the time instants these centers are achieved. 

Usage: 

Call CELL '(KSCAN,THM,OV, MOTION, VANG, VLIN,XKS, 

SWTCH , NCELL , XCTR , YCTR , TCTR , ELEV , AZIM , 

NCELLD) 


Input Parameters: 


KSCAN 


THM 
OV . 
MOTION 


VANG 

VEIN 

XKS 

SWTCH 

- NCELLD 


TYPE OF TRAJECTORY 

1 - SQUARE 

2 - HEXAGONAL 

3 - SPIRAL 

RADIUS OF THE UNCERTAINTY REGION (RADIANS) 
OVERLAP BETWEEN CELLS (%) 

DESIRED ANTENNA MOTION 

1 - CONSTANT ROTATIONAL SPEED 

2 - CONSTANT -SPEED ALONG TRAJECTORY 

3 - COMBINATION OF "1" AND "2" 

ANGULAR VELOCITY' (RADIANS/S). 

LINEAR VELOCITY ALONG TRAJECTORY (RADIANS/S) 
DISTANCE BETWEEN SUCCESSIVE CENTERS (RADIANS) 
SWITCHING ANGLE (IN RADIANS) TO GO FROM 
"MOTION=l" TO "MOTION=2" 

DIMENSION OF VECTORS XCTR, YCTR, TCTR, AZIM, AND 
ELEV IN THE MAIN PROGRAM 


Output Parameters: 


(XCTR, YCTR) 
TCTR. 

ELEV 

AZIM 


- VECTORS WITH THE CENTERS OF THE CELLS 

- VECTOR WITH THE TIME INSTANTS THE CENTER OF 
EACH CELL IS ACHIEVED 

-VECTOR WITH THEELEVATION OF THE CELLS CENTERS 

- VECTOR WITH THE AXIMUTE OF THE CELLS CENTERS 


Subroutines Required: 


QPOSTN,HPOSTN,SPOSTN 


6. Subroutine QPOSTN 

^ U 

Purpose: Determine the (K+1) — cell for the square type 

of trajectory. 

Usage: Call QPOSTN (K) 
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input Parameters: 

K - CELL NUMBER - 1 

ELI,AZI - ELEVATION AND AZIMUTH OF PREVIOUS CELL (COMMON) 

THB - RADIUS OF ANTENNA MAIN LOBE (COMMON 

XI, YI - CARTESIAN COORDINATES OF PREVIOUS CELL (COMMON) 

Output Parameters: 

ELO,AZO - ELEVATION AND AZIMUTH OF FOLLOWING CELL (COMMON). 

XO,YO - CARTESIAN COORDINATES OF FOLLOWING .CELL .(COMMON) 

Subroutines Required: None 

Subroutine HPOSTN 

Purpose: Determine the (K+l)~cen for the hexagonal type 

• ^ of trajectory. 

Usage: Call HPOSTN (K) 

Input Parameters: 

K - CELL NUMBER. - 1 

ELI,AZI - ELEVATION AND AZIMUTH OF PREVIOUS CELL (COMMON) 
THB - RADIUS OF ANTENNA MAIN LOBE (COMMON) 

Output Parameters: 

ELO,AZO - ELEVATION AND AZIMUTH OF FOLLOWING CELL (COMMON:) 

XO,YO - CARTESIAN COORDINATES OF FOLLOWING CELL (COMMON) 

Subroutines Required: None 

Subroutine SPOSTN 

■f* h 

Purpose: Determine the (K+l)“cell for the spiral type 

of trajectory. 

Usage: Call SPOSTN (K) ORIGINAL PAGE IS 

Input Parameters : POOR QUALIIY 

K - CELL NUMBER - 1 

ELI,AZI - ELEVATION AND AZIMUTH OF PREVIOUS CELL (COMMON) 
THB - RADIUS OF ANTENNA MAINLOBE (COMMON) 

XI, YI - CARTESIAN COORDINATES OF PREVIOUS CELL (COMMON) 


Output Parameters: 

ELO,AZO - ELEVATION AND AZIMUTH OF FOLLOWING CELL. (COMMON) 

XO,YO - CARTESIAN COORDINATES OF FOLLOWING CELL (COMMON) 
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Subroutines Required: ITR 
Subroutine ITR 

Purpose: Determine iteratively the azimuth of the next 

center, given the position of the previous one, 
for the spiral trajectory. 

Usage: Call ITR(A,B,AZO,GUSS,THB,D,PI ,K) 

Input Parameters: 


A,B - ELEVATION AND AZIMUTH OF PREVIOUS CENTER 

GUSS - INITIAL AZIMUTH GUESS 
THB RADIUS OF ANTENNA MAIN LOBE 

D - HALF DISTANCE BETWEEN SUCCESSIVE CENTERS 

PI _ - CONSTANT (3.141592653) 

K - CONTROL PARAMETER 

1 - INTERMEDIATE RESULT 
• 2 - FINAL RESULT 

Output Parameter: 

AZO - DESIRED AZIMUTH 

Subroutines Required: None 

. Subroutine TCENTR 


Purpose: Determine the' time instant. The cell (EL0,AZ0) 

is achieved following the normal path. 

Usage: Call TCENTR (KSCAN,M0TI0N,VANG,VLIN,D2,SWTCH) 


Input Parameters: 

KSCAN 

- TYPE OF TRAJECTORY 


1 - SQUARE 


2 - HEXAGONAL 


3 - SPIRAL 

MOTION 

- TYPE OF ANTENNA MOTION 


1 - CONSTANT ANGULAR VELOCITY 


2 -■ CONSTANT VELOCITY ALONG TRAJECOTYR 


3 - COMBINATION OF "1" AND '"2" 

VANG 

- ANGULAR VELOCITY (RADIANS/S) 

VLIN 

- LINEAR VELOCITY ALONG TRAJECTORY (RADIANS/S) 

D2 

- DISTANCE BETWEEN SUCCESSIVE CENTERS 

SWTCH 

- SWITCHING ANGLE (IN RADIANS) FROM "MOTION=l" 


"M0TI0N=2" 


TO 
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Output Parameter: 

TO • - INSTANT OF TIME A GIVEN CENTER IS ACHIEVED (COMMON) 
Subroutine Required: None 

11. Subroutine SLCELL 


Purpose: Determine number and position of sidelobe or nearby 

cells when the antenna is pointing to the origin 
(1 CELLED). 

Usage: Call SLCELL(KSL,THB,OVLAP,VANG,XSL,YSL,NSL,DTSL, 

NSLDIM) 


Input Parameters; 

KSL SELECTS SIDEL0BE(=1) OR NEARBY(=0) CELLS 

THB - BEAMWIDTH RADIUS (RADIANS) 

OVLAP - OVERLAP {%) 

VANG - ANGULAR VELOCITY (RADIANS/S) 

HSLDIM - DIMENSION OF VECTORS XSL AND YSL IN THE 

CALLING PROGRAM 


Output Parameters; 

NSL - NUMBER OF SIDELOBE CELLS 

(XSL,YXL) - VECTORS WITH THE CENTERS OF THESE SIDELOBE CELLS 

WHEN THE ANTENNA IS POINTING TO THE ORIGIN 

DTSL - TIME INTERVAL TO GO FROM ONE SIDELOBE CELL TO 

THE NEXT ONE 

Subroutines Required; None 


12. Subroutine TDRSS 

Purpose: Generate a TDRSS position in the uncertainty region 

according to a trauncated Gaussian probability density 
function 


Usage: 


P ( X , Y ) = KP*EXP ( - ( X**2+Y**2 ) / ( 2*SI GMAP**2 ) / ( 2 ( P I*S I GMAP**2 ) , 
X**2+Y**2.LE.THM**2, 

=0, ELSEWHERE, 

WITH 

KP= ( 1 - EXP ( -THM**2/ ( 2*S IGMAP**2 )))**(-!) 

Cal 1 TDRSS ( S0 ,SI GMAP , THM2 , XTDRSS , YTDRSS ) 


/-tTj-TnTNT&T, PAGE IS 
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Input Parameters: 

SO - VARIABLE USED IN THE RANDOM GAUSSIAN 

VARIABLE GENERATION 

SIGMAP - DESIRED STANDARD DEVIATION 

THM2 . - RADIUS OF UNCERTAINTY REGION SQUARED 

Output Parameters: 

(XTDRSS,YTDRSS)~ GENERATED TDRSS POSITION 
Subroutines Required: GAUSS 


. Subroutine GAUSS 


Purpose: Generate "N" (even) normally distributed random 

numbers “X" with mean "AVR" and standard deviation 
, "SIGMA". 

Usage: Call GAUSS(S0,N, AVR, SIGMA, X) 

Input Parameters: 

S0 - INPUT RANDOM NUMBER 

N - NUMBER OF SAMPLES 

AVR - DESIRED MEAN 

SIGMA - DESIRED STANDARD DEVIATION 

Output Parameters: 

X - VECTOR OF LENGTH "N" WITH THE GENERATED 

NORMALLY DISTRIBUTED NUMBERS 

Subroutines Required: UNIFOR 

Method: Box and Muller 

If 

XI, X2 - INDEPENDENT RANDOM VARIABLES UNIFORMLY DISTRIBUTED 

BETWEEN 0 AND 1 . 

THEN 

fl = SIGMA*SQRT(-2*AL06(X1))*C0S(2*PI*X2)+AVR, 
f2 = SIGMA*SQRT(-2*AL0G(X1 })*SIN(2*PI*X2)+AVR 
ARE NORMALLY DISTRIBUTED ACCORDING TO N(AVR, SIGMA) 
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SUMMARY 


This document presents the results pertaining lo the trade- 
off analysis and performance of the Ku-Band Shuttle antenna 
pointing and signal acquisition system. The study was performed 
assuming the existence of various antenna scanning trajectories 
and various signal acquisition algorithms. The square, hexagonal 
and spiral trajectories were investigated assuming the TORS 
postulated uncertainty region and a flexible statistical model 
for the location of the TORS within the uncertainty volume. 

The scanning trajectories, Shuttle/TDRS signal parameters and 
dynamics and three signal acquisition algorithms were integrated 
into a hardv^are simulation discussed herein and documented in 
detail in Volume IV- The hardware simulation is quite flexible 
in that it allows one to evaluate signal acquisition -performance . 
for. an arbitrary (programmable) antenna pattern, a’ large range' 
of C/Ng‘s, various TDRS/Shuttl e. a priori' uncertainty distributions 
and three distinct signal search algorithms. 

Based upon the data made available during this contract 
period, certain Ku-Band forward link signal threshold characteristics 
were studied. The antenna pointing and acquisition threshold is found 
to be less than 60 dB-Hz vn’th an acquisition time dependent on the 
antenna scan procedure and acquisition algorithm implemented. 

Various techniques are discussed in this report and a computer 
program is presented from which these can be evaluated. It 
appears that there will be no problem in meeting t!ie system 
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spec even in light of the uncertainty associated v/ith the 

possibilities of antenna sidelobe acquisition. The details ffe 

provided herein. Assuming a single channel monopulse system and a 

30 MHz IF bandwidth, a monopulse tracking loop bandwidth of one Hz, 

the monopulse tracking jitter is 0.11 degrees at C/Nq = 54 dB-Hz. 

The Costas loop arm filter bandwidths can be chosen to be ten 

times the data rate such that no false-lock problem occurs 

during carrier acquisition. Assuming loop bandwidth of 

= 3 kHz, the loop jitter is 10 degrees at C/Nq = 60.4 dB-Hz. 

Acquisition can -be accomplished in less than 10 seconds. 

For the purpose of antenna scanning analysis, a reference 

coordinate system, whose z-axis is in line with the center 

axis of a specified uncertainty cone of the TORS position, is 

chosen for the relative relation between the TORS and the Shuttle 

antenna. The scan path of the Shuttle antenna can be projected 

onto the (x,y) plane of the coordinate system; v/hile the uncertainty 

■ cone of tha TORS position can be described by a circle in the 

(x,y) plane. The uncertainty in the position of the TDRS is 

modeled by a truncated Gaussian probability density p(x,y) with 

2 2 2 . ' 

uncertainty parameters a = a = a . By changing the value 

X j 

of this variance parameter, the model is sufficiently general 
to include a uniform distribution of TDRS position uncertainty 
to one v/hich specifies the position with probability one 

(Op=0). 

In the study of antenna scanning three types of trajectories 
are proposed and evaluated in terms of average scan time and the 
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structure of coverage over the uncertainty cone of the TORS 
position. The analysis technique used for finding average 
scan time is discussed and an illustration of the technique 
for a specific type of scanning trajectories and motion of the 
Shuttle antenna is given. The results show that the spiral tra- 
jectories is in general better than the other two trajectories* 
especially for a constant velocity along a trajectory, since the 
path of a spiral trajectory (from the center of uncertainty cone 
to its edge) can be shorter than the other two. 

From the Ku-band system specifications, it has been 
recognized that the variation of the received signal level at 
the Shuttle (due to TORS EIRP path loss and antenna pointing loss 
variations) varies as much as 23 dB. If the antenna sidelobes are 
not sufficiently suppressed, a potential problem called sidelobe 
acquisition may cause the degradation in system performance. 
Therefore in this tradeoff analysis study, three different 
acquisition strategies are -proposed and evaluated. The first 
acquisition strategy may be used for the case where the sidelobe 
acquisition does not impose a problem to the system. The other 
two strategies are primarily designed to avoid the sidelobe 
acquisition (especially the first sidelobe). 

In this study, wc assume that the first sidelobe of the 
Shuttle antenna is 17.5 dB suppression from the main lobe of the 
antenna so that v/e can study the effectiveness of the acquisition 
algorithm proposed. Surely, one should note that the sidelobe 
acquisition can be overcome by tapering down the sidelobes of 
the Shuttle antenna below 23 dB. However, it has been 
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found that with one particular acquisition algorithm studied, 
a strong sideTobe of the antenna can speed spatial acquisition 
significantly. 

The scan schemes and acquisition algorithms are all integrated 
by a computer simulation program. This software package, discussed 
herein, provides a useful tool for predicting the performance 
of the Ku-band antenna pointing system. Various options are 
available for users to do tradeoff studies on system parameters 
in designing an acquisition system. The software package has 
been tested on UNIVAC 1100 Series computers and verified with 
analytical results. 
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